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ABSTRACT 

A comparative study of spectra of 21 Type la supernovae (SNe la) obtained 
about one week before maximum light, and 8 spectra obtained 11 or more days 
before maximum, is presented. To a large extent the premaximum spectra ex- 
hibit the defining characteristics of the four groups defined in Paper II (core- 
normal, broad-line, cool, and shallow-silicon). Comparisons with SYNOW syn- 
thetic spectra show that all strong features and most weak ones can be accounted 
for in a plausible way. The issues of detached high-velocity features, the possible 
ubiquity of carbon clumps, the maximum detectable ejecta velocities, and the 
possibility of blueshifted emission-line peaks are discussed. 

Subject headings: supernovae: general 

1. INTRODUCTION 

This is the third in a series of papers on a comparative direct analysis of the optical 
spectra of Type la supernovae (SNe la). Paper I (Branch et al. 2005) was concerned 
with a time series of spectra, from 12 days before to 115 days after the time of maximum 
light, of the well-observed, spectroscopically normal, Type la SN 1994D. Paper II (Branch 
et al. 2006) concentrated on near-maximum-light spectra of 24 SNe la. In this paper we 
focus on pre-maximum-light spectra, which form in the outermost layers of the ejected 
matter and provide clues to the nature of the explosions that cannot be obtained from later 



1 Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma, Norman, OK 73019; 
e-mail: branch@nhn.ou.edu 

2 Division of Network Business, PionetSoft Corporation, Shibuya-ku, Tokyo 151-0072, Japan 

3 Department of Astronomy, Indiana University, Bloomington, IN 47405 

4 Department of Astronomy, Whitman College, Walla Walla, WA 99362 



- 2 - 



spectra that form deeper in the ejecta. Our emphasis is on line identifications, which are 
needed in the quest to develop good hydrodynamical explosion models for SNe la. When 
detailed spectrum calculations are carried out for explosion models (e.g., Baron et al. 2006; 
Howell et al. 2001), discrepancies between synthetic and observed spectra invariably appear. 
Knowing the identities of the observed features provides clues to how the explosion models 
should be modified in order to achieve better agreement with observation. 

In Paper II we divided the maximum-light spectra into four groups: core-normal; 
broad-line; cool; and shallow-silicon (denoted CN, BL, CL, and SS, respectively, in the 
figures and tables of this paper). The group assignments were made on the basis of mea- 
surements of the (pseudo) equivalent widths of absorption features near 5750 A and 6100 A, 
as well as on the appearance (depth, width, and shape) of the 6100 A absorption, which is 
produced by the Si II A6355 transition. Although we framed the presentation and discussion 
in terms of the four groups, in the end we concluded that for the most part the spectra 
appeared to have a continuous distribution of properties, rather than consisting of discrete 
subtypes. Therefore, exactly how the group boundaries were drawn was not important. 

Figure 1, a plot of the equivalent width H^(5750) of the absorption feature near 5750 A 
against the equivalent width 1^(6100) of the absorption feature near 6100 A m maximum- 
light spectra, is an update of the corresponding figure in Paper II, with the following addi- 
tional events: the core-normal SN 2003cg (Elias-Rosa et al. 2006), SN 2003du (Gerardy et 
al. 2004), and SN 2004S (Krisciunas et al. 2006); the broad-line SN 1983G (McCall et al. 
1984) and SN 2002er (Kotak et al. 2006); and the shallow-silicon SN 1999aa (Garavini 
et al. 2004), SN 1999ac (Garavini et al. 2005), SN 2000E (Valentini et al. 2003), SN 2005cg 
(Quimby et al. 2006), and SN 2003fg [the proposed super-Chandrasekhar SN la, also known 
as SNLS-03D3bb; Howell et al. (2006); Jeffery, Branch, & Baron (2007a); but see also Hille- 
brandt, Sim, & Ropke (2007)]. SN 1999ac is a borderline case: in Figure 1 it falls between 
two core-normals but on the basis of the shape, depth, and width of the 6100 A absorption, 
we include it with the shallow-silicons. 

To a large extent the premaximum spectra turn out to exhibit the defining characteristics 
of the four maximum-light groups of Paper II, so we use the same terminology in this paper 
(with all group assignments based on maximum-light, not premaximum, spectra). Again, we 
are not suggesting that SNe la actually consist of discrete subtypes; instead, we are thinking 
in terms of three inhomogeneous groups radiating outward (in Figure 1) from the highly 
homogeneous core normals. 

In Figure 1, three lines of constant H^(5750)/TV(6100) ratio are shown. This ratio is 
similar to although not the same as the well known R(Si II) ratio (Nugent et al. 1995), 
which is based on the depths of the same two absorption features. Considerable diversity 



- 3- 



at a given W (57 50) /W (6100) ratio is evident, especially at the low value of 0.05, e.g., the 
SN 2005hk, SN 1999ee, and SN 2002bf have similar W (5750) /W (6100) ratios but otherwise 
thei spectra are very different. 

In this paper, as in Papers I and II, we confine our attention to optical spectra, from 
the Ca II H&K feature in the blue (~ 3700 A) to the Ca II infrared triplet (Ca II IR3) in 
the red (~ 9000 A). All spectra have been corrected for the redshifts of the parent galaxies. 
The original spectra have a range of slopes, owing to intrinsic differences, differences in the 
amount of reddening by interstellar dust, and observational error. Because we are interested 
only in the spectral features, and not in the underlying continuum slopes, all spectra in this 
paper, both observed and synthetic, have been flattened according to the local normalization 
technique of Jeffery et al. (2007b). The flattening facilitates comparison of the spectral 
features. Mild smoothing also has been applied to some of the spectra. 

We examined all of the SN la premaximum spectra available to us and selected two 
samples: a "one- week premax" sample, consisting of one spectrum each of 21 SNe la observed 
between day —5 and day —9 with respect to the time of maximum brightness in the B band, 
and an "early" sample, consisting of one spectrum each of eight SNe la observed on or earlier 
than day —11. (All SNe of the early sample are also in the one-week premax sample.) The 
SNe la and the epochs of the selected spectra are listed in Table 1. 

Continuing our attempt to provide an internally consistent quantification of SN la spec- 
tra, we have used the SYNOW code to fit all spectra of Table 1. Descriptions of SYNOW 
and its use can be found in Branch et al. (2003) and Paper fj]. Compared to other supernova 
synthetic-spectrum codes, SYNOW is simple, but when working in an empirical spirit from 
the data, as opposed to evaluating an explosion model, SYNOW offers the advantage that 
the user has direct control over line optical depths as function of radial coordinate, while 
other codes input composition and density structures and then compute the temperature 
structure and the line optical depths. Thus SYNOW is the most useful code for making line 
identifications (pending eventual confirmation by means of detailed calculations for explo- 
sion models whose spectra and light curves match the observations). For the present work, 
we adopt precepts for fitting spectra that are much like those of Paper II, e.g., the default 
excitation temperature is 7000 K for the cool SN 1986G and SN 1999by and 10,000 K for 
all others. 



1 For consistency with Papers I and II, in this paper we do not use SYNOW version 2.0, described and 
used by Parrent et al. (2007). 
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2. CORE-NORMALS 

Figure 2 shows the spectra of the seven core-normals of the one-week premax sample. 
The spectrum of SN 2001el differs from the others; as shown by Mattila et al. (2005), 
it resembles the day —14 spectrum of SN 1990N, so it will be discussed below, with the 
early sample. The other spectra show considerable homogeneity, although the diversity is 
somewhat greater than in the maximum-light core-normal sample. The 6100 A absorption 
is deeper in SN 1998bu and SN 2003cg than in the others, and the broad absorption near 
4300 A is weaker in SN 1994D and SN 1998aq than in the others. SN 2003du has strong 
HV§ Ca II absorptions near 3700 A and 8000 A (although not nearly as strong as in the 
exceptional SN 200 lei). SN 1990N also has a strong absorption near 3700 A. Diversity m 
HV Ca II was present also in maximum-light spectra of core-normals. 

The SYNOW fitting parameters for core normals of the one-week premax sample are in 
Table 2. An example fit, to the day —7 spectrum of SN 2003du, is shown in Figure 3. The 
fit is good and most of the line identifications seem clear. The most obvious discrepancy, the 
failure of the synthetic spectrum to match the observed absorption near 3980 A, is primarily 
caused by Si II A4130, so we can not remove this discrepancy without doing harm to the fit 
elsewhere. 

The synthetic absorption produced by Si III A4560 is a bit too blue to match the observed 
absorption at 4400 A (which is present in most observed spectra of this paper). A better 
fit could be achieved by using C III A4649 instead of Si III. But in a synthetic spectrum 
that contains only Si III lines the A4560 absorption does occur at the right place to account 
for the 4400 A absorption, and Si III A5743 also accounts for a weak observed absorption 
near 5560 A. In most premaximum spectra there is little or no evidence for C II features, 
so it would be surprising to see ubiquitous evidence for the high-excitation (29.6 eV) C III 
line. Therefore, even though neither Si III feature fits well in the full synthetic spectrum 
(the A5743 absorption does not even appear as a distinct feature), we assume that these two 
Si III lines usually are responsible for the corresponding observed features. 

The synthetic absorption produced by Si II A5972 does not account well for the observed 
absorption at 5750 A. This and the other two discrepancies mentioned above are generic to 
our fits. The fit could be improved by invoking Fe III lines with a higher excitation potential, 
at or approaching 15,000 K, so that lines near 6000 A contribute to the 5750 A absorption 
[see Branch et al. (2003) on SN 1998aq]. But the problem with the 5750 A absorption 



2 As in Papers I and II, PV refers to features that form at or near the photospheric velocity and HV refers 
to high-velocity detached features. 
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occurs in some spectra in which we see no evidence for the stronger Fe III lines in the blue, 
so we refrain from using high-excitation Fe III to resolve this discrepancy. Another way to 
improve the fit, in some cases, would be to use a higher Si II optical depth and impose a 
maximum velocity on Si II (see §3 and §4). 

In our fits to the core normals of the one-week premax sample, the relative weakness 
of the absorption near 4300 A in SN 1994D and SN 1998aq is accounted for by using lower 
Mg II and Fe III optical depths (see Table 2). In the same two SNe we have used C II lines, 
to allow A6580 to improve the fit near the peak of the Si II A6355 feature. 

Figure 4 shows the spectra of the three core-normals (and the two broad-lines, to be 
discussed in §3) of the early sample. The diversity among the core-normals of the early sam- 
ple is much greater than in the one-week premax sample, e.g., the three 6100 A absorptions 
have distinctly different shapes, and in SN 1994D the absorption trough from 4600 to 5000 A 
is deeper than in SN 1990N and SN 2003du. The early spectrum of SN 1994D is an example 
of group crossover: the spectrum is more like that of the broad-line SN 2002bo than like the 
other core-normals. (On the other hand, the spectrum of the broad-line SN 2002er could 
pass for core-normal.) 

The SYNOW fitting parameters for the three core-normals (and the two broad-lines) 
of the early sample are in Table 3. The deeper 4600 A to 5000 A absorption trough in 
SN 1994D is accounted for primarily by using stronger HV Fe II. 

How to fit the various shapes of the 6100 A absorption is an interesting issue. For 
SN 1994D we fit it in the conventional way, with a simple PV Si II component (a fit is 
shown in Paper I). For SN 2003du, the feature appears to consist of two marginally resolved 
components (Stanishev et al. 2007), so we use both PV and HV Si II. The fit is shown in 
Figure 5. (We did not need to use HV Si II for the one-week premax spectrum of SN 2003du.) 

The unusual, nearly flat-bottomed 6100 A absorption of the early spectrum of SN 1990N 
has been discussed by several authors (Jeffery et al. 1992; Fisher et al. 1997; Mazzali 2001). 
Mattila et al. (2005) showed that the one-week premax spectrum of SN 2001el resembled 
the early spectrum of SN 1990N (see their Figure 1). As can be seen by comparing our 
Figures 2 and 4, the HV absorption of the Ca II infrared triplet is much deeper in SN 2001el 
(requiring a covering factor approaching unity), but otherwise the two spectra, including the 
6100 A absorption, are quite similar. Mattila et al. attributed the flat-bottomed shape of 
the 6100 A absorption to line formation in a shell that is geometrically thin compared to the 
size of the photosphere, but this explanation is not viable. Line formation in a thin shell that 
is not detached from the photosphere does produce a (shallow) flat-bottomed absorption [see 
Figure 8 of Jeffery & Branch (1990)], but the blue edge of the absorption is blueshifted only 
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by an amount that corresponds to the photospheric velocity, which is not the case in the 
spectra of SN 1990N and SN 2001el. We can nicely fit the flat-bottomed absorptions with a 
single Si II component, using a high value of v e and imposing a maximum velocity. However, 
since the blue edge of the absorption corresponds to a high velocity comparable to that at 
which the HV Ca II and HV Fe II are forming (~ 20, 000 km s _1 ), it seems more reasonable 
to assume that the the feature is is an unresolved blend of PV and HV components, and fit 
the absorption accordingly (see Table 2 for SN 2001el and Table 3 for SN 1990N). Support for 
this interpretation is provided by a previously unpublished spectrum of SN 1990N obtained 
at day —12 in bad weather by S. Benetti and M. Turatto (Figure 6). The absorption feature 
resembles the marginally-resolved two-component absorptions of SN 2003du (Figure 5) and 
SN 2005cg (Figure 13). [see also Garavini et al. 2007 on 05cf] 



3. BROAD-LINES 

Figure 7 shows the spectrum of the core-normal SN 2003du (for comparison) and the 
spectra of the four broad-lines of the one-week premax sample. In maximum-light spec- 
tra (Paper II) broad-lines have the same spectral features and line identifications as core- 
normals, but the equivalent width of the 6100 A absorption is larger and the features are 
generally broader. The four broad-lines of Figure 7 also have these characteristics, although 
as in Paper II, the borderline case SN 1992A differs from core-normal mainly in the depth 
and equivalent width of the 6100 A absorption. SN 2002er is much like SN 1992 A except 
for a weaker absorption near 7450 A. As at maximum light, SN 1984A is the most extreme 
broad-line in the sample. 

The SYNOW fitting parameters for broad-lines of the one-week premax sample are 
listed in Table 4. As at maximum light (Paper II), the one-week premaximum spectra of 
the broad-lines can be fitted with the same ions as used for the core-normals, but with 
higher values of v e and/or rl Figure 8 shows our fit to the day -7 spectrum of the extreme 
broad-line SN 1984 A. For Si II a high value of v e = 4000 km s _1 is used and a maximum 
velocity of 27,000 km s _1 is imposed. 

The spectra of the two broad-lines of the early sample are shown in Figure 4, and fitting 
parameters are in Table 3. Fe III is not used for the early spectra of the broad-lines, but its 
features could be hidden by the strong HV Fe II lines. Figure 9 shows our fit to the day —14 



3 When inspecting fitting parameters one should keep in mind the partial degeneracy of v e and r. For 
example, for SN 2002bo the synthetic absorption produced by Si II A6355 with r = 7 and v e = 2000 km s _1 
is roughly similar to the corresponding absorption in SN 1992A with r = 45 and v e = 1000 km s -1 . 
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spectrum of SN 2002bo. In an attempt to account for the 5750 A absorption with Si II A5972 
we have used a large optical depth of 30 for Si II and imposed a maximum velocity of 25,000 
km s~ x . A maximum velocity of 22,000 km s -1 has been imposed on S II. 

4. COOLS 

Figure 10 shows the spectrum of the core-normal SN 2003du and the spectra of the 
three cools of the one-week premax sample. For reasons given in Paper II, the borderline 
SN 1989B was included in the cool group even though it did not show the blue trough from 
about 4000 A to 4300 A that is characteristic of the other cools. Here again SN 1989B has 
a higher H^(6100) value than core-normals and resembles SN 1986G at wavelengths longer 
than about 4400 A, while lacking the blue trough. SN 1989B remains a borderline case. 

The SYNOW fitting parameters for cools of the one-week premax sample are listed 
in Table 5. At maximum light (Paper II) our strategy when attempting to simultaneously 
account for the 5750 A and 6100 A absorptions was to use a high optical depth and impose 
a maximum velocity on Si II (as we have done in this paper for some of the broad-lines; see 
§3). Figure 11 shows our fit to the day —5 spectrum of SN 1999by, where we have used the 
same strategy with only partial success; in the synthetic spectrum, the 5750 A absorption 
is too weak even though the 6100 A is too strong. We regard the line identifications in 
Figure 11 to be definite, except for Sc II, which helps in one place but cannot be considered 
definite. 

There are no cools in the early sample. Only a small number of cools have been observed 
at all, and because they have short rise times and are dim relative to other SNe la, no spectra 
have been obtained at times as early as day —11. 

5. SHALLOW-SILICONS 

Figure 12 shows the spectra of the core-normal SN 2003du and the seven shallow-silicons 
of the one-week premax sample. These were classified as shallow-silicon on the basis of their 
maximum-light spectra but they also have shallow-silicon at one-week premax. (SN 1999ac 
is borderline in this respect but it is more like the shallow-silicons than the core normals 
in other respects, such as in having weak S II absorptions.) There is plenty of diversity, 
however, extending to the extreme cases of SN 1991T and its near twin SN 1997br which at 
a glance do not necessarily have Si II at all (although we do include it in our fits), and to 
SN 2005hk. SN 2005hk is in our one-week premax sample of shallow-silicons (with fitting 
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parameters in Table 6 the same as used in Chornock et al. 2006), and SN 2002cx was in our 
maximum-light sample of shallow-silicons. But as discussed in Paper II and elsewhere (Li 
ct al. 2003; Branch et al. 2004b; Jha et al. 2006; Chornock et al. 2006; Stanishev et al. 
2007; Phillips et al. 2007), in spite of certain spectroscopic similarities to SN 1991T-likes, 
the SN 2002cx-likes appear to be physically distinct from other SNe la. 

The fitting parameters for the shallow-silicons of the one-week premax sample are in 
Table 6. The shape of the 6100 A absorption in SN 2005cg is especially interesting, because 
it appears to consist of marginally resolved PV and HV components of Si II. Our fit is 
shown in Figure 13. Note that in Figure 12 the 6100 A absorption of SN 1999ee bears some 
resemblance to that of SN 2005cg. As with the one-week premax spectrum of SN 2001el 
and the early spectrum of SN 1990N, the 6100 A absorption of SN 1999ee can be fit with 
a single broad PV component or an unresolved combination of a PV and an HV feature. 
For the same reasoning as in the discussion of SN 1990N and SN 2001el and the results of 
Garavini et al. (2007) and Stanishev et al. (2007), we use the two components. 

Figure 14 shows the spectra of the core-normal SN 2003du and the three shallow- 
silicons of the early sample. Again SN 1999aa and SN 199 IT do have shallow silicon, and 
again SN 1999ac is borderline. Figure 15 shows our fit to the day —12 spectrum of SN 1999aa. 
Like Garavini et al. (2004), we have resorted to C III, but the identification is not definite. 
When fitting the HV Ca II infrared triplet absorption, the HV Ca II H&K absorption comes 
out too deep; this happens in some of our other fits, although less than here. We believe 
the absorption feature near 5500 A to be Si III A5743, but as usual it does not appear as a 
distinct feature in the synthetic spectrum. 

6. DISCUSSION 

The previous sections have been organized in terms of the four groups of Paper II, with 
all group assignments based on the appearance of near-maximum-light spectra. To a large 
extent the one-week premax spectra of this paper retain the defining characteristics of the 
four groups, i.e., group assignments based on one-week premax spectra would be much the 
same. To a lesser extent, the same is true of the early spectra of this paper, e.g., almost 
all shallow-silicons are distinct from core-normals (Figure 14). But there is some crossover 
between groups in the early sample, e.g., Figure 4 shows that the early spectra of the core- 
normal SN 1994D and the broad-line SN 2002er are not very different. In any case, as in 
Paper II we see little or no evidence that SNe la actually break up into distinct groups, 
except for the SN 2002cx-likes (and possibly the cools, but see Pastorello et al. 2007 on 
SN 2004eo for evidence of a connection between core normals and cools). 
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The SYNOW fits to the premaximum spectra of this paper generally are of the same 
quality as the fits to the maximum-light spectra of Paper II. All strong features and most 
weak ones can be accounted for in a plausible way. Most line identifications seem clear, 
although as discussed in previous sections there are some ambiguities, and some of the 
weak features remain unidentified. For the most part our line identifications in premaximum 
spectra are consistent with those suggested previously on the basis of SYNOW (or SYNOW- 
like) synthetic spectra, as well as synthetic spectra calculated with the Mazzali-Lucy Monte 
Carlo code. Such studies have been carried out for SN 1990N (Jeffery et al. 1992; Mazzali 
et al. 1993; Mazzali 2001), SN 1991T (Ruiz-Lapuente et al. 1992; Jeffery et al. 1992; 
Mazzali, Danziger, & Turatto 1995; Fisher et al. 1999), SN 1997br (Hatano et al. 2002), 
SN 1998aq (Branch et al. 2003), SN 1999aa (Garavini et al. 2004), SN 1999ac (Garavini 
et al. 2005), SN 1999ee (Mazzali et al. 2005a), SN 2002bo (Benetti et al. 2004; Stehle et al. 
2005), SN 2003cg (Elias-Rosa et al. 2006), and SN 2005hk (Chornock et al. 2006). 

One of the notable differences from line identifications in some of the previous papers is 
that neither in this paper nor in Paper II do we use PV Fe II or Co II (except in SN 2005hk). 
In Paper I on SN 1994D we first used weak PV Fe II at day +4, and at later epochs it became 
much stronger. This appears to be consistent with the argument of Kasen & Woosley (2007) 
that the development of extensive line blocking by PV Fe II (and Co II) lines soon after 
maximum light is the key to understanding the SN la width-luminosity relation (Phillips 
et al. 1999). We will return to this issue in Paper IV, on post-maximum SN la spectra. 

We do use PV Fe III for the one-week premax sample, except in the cools. It is commonly 
said, correctly, that SN 1991T-likes are characterized by conspicuous Fe III features. It is 
not true, though, that the Fe III features are particularly strong in SN 1991T-likes. At 
one-week premax we use Fe III optical depths for the shallow-silicons that are comparable 
to or lower than for the core-normals and the broad-lines. The same was true at maximum 
light (although we failed to mention it in Paper II). At maximum and one-week premax, 
Fe III features are conspicuous in shallow-silicons only because there is less competition from 
other features. 

It is interesting to consider some of the implications of the fitting parameters, especially 
the line optical depths, in the context of the local-thermodynamic-equilbrium (LTE) calcu- 
lations of Hatano et al. (1999a). For example, Figure 5 of Hatano et al. shows that in LTE 
near 7500 K the Si II optical depth can exceed unity in a composition in which hydrogen 
and helium have burned to carbon and oxygen while the mass fraction of other elements, 
including silicon, is only solar. This suggests that the presence of Si II in the spectrum 
does not necessarily require synthesized silicon. However, for this composition the optical 
depth of Fe II is predicted to be greater than that of Si II, which according to us is not 
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the case in premaximum or near-maximum spectra. Figure 6 of Hatano et al. shows that 
for a carbon-burned composition, which does include synthesized silicon, the Si II optical 
depth exceeds that of Fe II by a large factor. Therefore, since we do not use PV Fe II in 
premaximum and maximum spectra, we expect that wherever PV Si II is detected, it is pro- 
duced by synthesized, rather than primordial, silicon. Another interesting example is that 
for the core-normals of the one-week premax sample, the Si II optical depths range from 
four to eight (and for the core-normals of the maximum-light sample it ranged from six to 
13). Figure 6 of Hatano et al., for the carbon-burned composition, shows that a temperature 
difference of only 800 K is sufficient to change the Si II optical depth by a factor of two. 
This temperature sensitivity makes the high degree of homogeneity among the core-normals 
all the more remarkable. 

The HV features are intriguing. Following the discovery of HV Ca II and HV Fe II in 
SN 1994D (Hatano et al. 1999b), HV Ca II has been recognized by many authors: see Li 
et al. (2001), Thomas et al. (2004), and Branch et al. (2004a) on SN 2000cx; Wang et al. 
(2003), Kasen et al. (2003), and Mattila et al. (2005) on SN 2001el; Mazzali et al. (2005a) on 
SN 1999ee; Gerardy et al. (2004) and Stanishev et al. (2007) on SN 2003du; Quimby et al. 
(2006) on SN 2005cg; and Garavini et al. (2007) on SN 2005cf. In premaximum spectra of 
SNe la, HV Ca II is ubiquitous (Mazzali 2005b; and this paper), and in Paper II and we 
regularly invoked HV Ca II even at maximum light. The identification of HV Fe II, which 
usually appears at modest strength in a crowded spectral region, is perhaps not definite, but 
we have invoked it regularly both in this paper and in Paper II at maximum light. One thing 
that is clear, from flux spectra (Thomas et al. 2004; Tanaka et al. 2006) and especially from 
polarization measurements (Wang et al. 2003; Kasen et al. 2003), is that in at least some 
cases the HV Ca II features form in asymmetrical structures. The origin of the HV features 
is unclear. One possibility is that for one reason or another they form naturally in the ejecta 
(Hatano et al. 1999b). An alternative, favored by Gerardy et al. (2004) and Quimby et al. 
(2006), is that they form in a shell of restricted velocity interval and elevated density that 
is produced by interaction between the ejecta and circumstellar matter (an accretion disk, 
a filled Roche lobe, or a common envelope). Our identification of HV Fe II, if correct, may 
hint that the features can form naturally in the ejecta, because more often than not we 
detached HV Fe II at a lower velocity than HV Ca II. In the spectra of the early sample, it 
is unlikely that the HV Fe II features are produced by iron resulting from decay the decay of 
56 Ni and 56 Co, because lines of HV Ni II and HV Co II also would be expected (see Figure 9 
of Hatano et al. 1999a), but it is possible that the HV Fe II features are produced by iron 
that has been freshly synthetisized along with 56 Ni. Another hint that the HV features may 
arise naturally in the ejecta is that in some cases we see HV Si II at about the same velocity 
as HV Ca II and HV Fe II. In the circumstellar interaction scenario, the blue wing of the 
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Si II absorption is predicted to be truncated rather than extended (Gerardy et al. 2004; 
Quimby et al. 2006). In any case, it is interesting that that HV features, of Ca II, Fe II, and 
Si II typically are detached not far from 20,000 km s -1 . The density and/or composition 
structure of most SNe la seems to have something special around 20,000 km s^ 1 , probably 
including a density increase (Mazzali et al. 2005a; Mazzali 2005b). 

The day —14 spectrum of the broad-line SN 2002bo (Figure 8 and Table 3) raises an 
issue about the definition of HV and PV features. We did not use HV Ca II or HV Si II, and 
what we called HV Fe II was only mildly detached, at 22,000 km s -1 , from the photosphere 
at 20,000 km s -1 . Thus the absence of HV Ca II or HV Si II is a matter of definition: Ca II 
and Si II are forming at 20,000 km s _1 , as in other SNe la, but they are not referred to as 
HV simply because of the high velocity at the photosphere. 

The issue of carbon also is interesting for constraining explosion models. As discussed 
in §2, a C III line could be used to improve our fits for core-normals but we suspect that it 
would be spurious. We have used C III to improve the fits for several shallow-silicons but in 
no case is the identification convincing. Marion et al (2006) have established that C I lines 
are not detected [except in cools: Howell et al. (2001)], but their conclusion that carbon is 
much less abundant than oxygen depends strongly on their assumption of LTE at 5000 K. A 
temperature of 5000 K seems too low to account for the strengths of other lines (e.g., Fe III) 
and at higher temperatures the upper limit on the carbon-to-oxygen ratio would be much 
lower (see Figure 5 of Hatano et al. 1999a). 

Lines of C II offer the best prospects for detecting the presence of carbon. Evidence 
has been presented for the strongest optical C II line, A6580, in SN 1998aq (Branch et al. 
2003); SN 1990N (Jeffery et al. 1992; Mazzali 2001); SN 1999ac (Garavini et al. 2005); and 
SN 1994D and SN 1996X (Branch et al. 2003 and Paper I). Thomas et al. (2007) presented 
the most convincing evidence, in SN 2006D. A A6580 absorption was clearly present at 
day —7 and day —5, and it weakened as maximum light was approached. In this paper we 
have invoked C II in these events and only one other: the early spectrum of SN 2003du, but 
in this case C II A6580 is being used mainly just to beat down the Si II A6355 emission peak 
and the identification is not convincing. Thus, even at premaximum, C II is elusive, or, as 
Thomas et al. put it, sporadic, perhaps due to its tendency to appear in clumps. If the 
clumps are few in number and none happen to be in front of the photosphere, no absorption 
features will be seen. Even if a clump is in front of the photosphere, if it has sufficiently 
high line-of-sight velocity (about 15,000 km s _1 or more) the A6580 absorption may be 
blueshifted into the Si II absorption and not be detectable. Therefore, although C II is seen 
in only a fraction of premaximum SN la spectra, the ubiquitous presence of carbon clumps 
in the outer layers of SNe la is not excluded. Unburned carbon at the observed velocities is 
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not characteristic of published delayed-detonation models. 

Another issue to be addressed on the basis of premaximum spectra is that of the max- 
imum detectable ejection velocities. From Figures 2, 4, 7, 10, 12, and 14 we can see that 
Ca II absorption is generally detectable up to about 25,000 km s _1 and sometimes up to 
about 30,000 km s^ 1 . The highest maximum velocity that has been imposed on Ca II is 
34,000 km s" 1 for SN 2001el. Detectable optical depth in Ca II can easily be produced with 
a solar mass fraction of Ca II, so the maximum velocities of synthesized Ca II may be lower. 
From the same figures we can see that Si II absorption is generally detectable up to 15,000 
km s _1 and sometimes above 20,000 km s -1 . The highest maximum velocity that has been 
imposed on Si II is 27,000 km s _1 for the one-week premax spectrum of SN 1984A. We have 
argued above that this absorption is produced by synthesized silicon. This is consistent with 
the 11,900 km s -1 that Mazzali, Benetti, & Hillebrandt (2007) adopted for a lower limit 
to the top of the silicon-rich layer, but the silicon-rich layer seems to generally extend to 
much higher velocities. Such high silicon velocities are not characteristic of pure deflagration 
models. The presence of carbon at intermediate velocities and silicon at high velocities could 
be consistent with the pulsating reverse detonation models of Bravo & Garcia-Senz (2006). 

Another velocity-related issue is the possibility that some of the P-Cygni emission 
components are blueshifted, as they sometimes are in SNe II (Chugai 1988; Dessart & Hiller 
2005). From measurements of the wavelengths of flux maxima Blondin et al. (2006) inferred 
emission-line blueshifts of up to 8000 km s" 1 in premaximum SN la spectra. The assump- 
tions on which SYNOW is based (sharp photosphere and resonant scattering) do not allow 
blueshifts of unblended emission-line peaks, yet our fits do not in general have a problem 
reproducing the wavelengths of the observed flux maxima. In SYNOW synthetic spectra the 
flux maxima are affected by line blending; sometimes they are a bit too red, consistent with 
an emission-line blueshift, but equally often they are a bit too blue. In SNe la, line blending 
is a serious problem for establishing that emission-line peaks are blueshifted. 

The premaximum spectra of SNe la exhibit a rich diversity, which even within each 
of our four groups probably is multidimensional. The outer layers of SN la ejecta differ in 
composition, density, and temperature structure, and they have asymmetries - probably both 
global shape asymmetries as well as composition clumping. Understanding the observational 
diversity in terms of the physical differences is a challenging task that will require more 
observations and modeling. Flux spectra have much to tell us, but in view of the asymmetries, 
polarization spectra are a high priority. 

We are grateful to Stefano Benetti and Massimo Turatto for providing the previously 
unpublished spectrum of SN 1990N, and to all observers who have provided their published 
spectra. This work has been supported by NSF grants AST-0204771 and AST-0506028, and 
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Fig. 1. — An update of Figure 1 of Paper II: W(5750) is plotted against W^(6100) for 
spectra obtained within three days of maximum light. Core-normal SNe la are shown as 
circles, broad-line SNe la as squares, cool SNe la as diamonds, and shallow-silicon SNe la 
as triangles. The labels of the dashed lines refer to the W {5750) /W '(6100) ratio. 
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Fig. 2. — Spectra of seven core-normals of the one-week premax sample. The spectra have 
been flattened as described in the text. Vertical displacements are arbitrary and narrow 
absorptions near 7600 A and 6900 A are telluric. Vertical dashed lines refer to Ca II A3945 
and A8579 blueshifted by 25,000 km s" 1 and Si II A6355 blueshifted by 15,000 km s _1 . 
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Fig. 3. — The day —7 spectrum of the core-normal SN 2003du (solid line), from Anupama 
et al. (2005), is compared with a synthetic spectrum (dashed line). 
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Fig. 4. — Spectra of three core-normals and two broad-lines of the early sample. Vertical 
dashed lines refer to Ca II A3945 blueshifted by 30,000 km s -1 , Si II A6355 blueshifted by 
20,000 km s" 1 , and Ca II A8579 by 25,000 km s" 1 . 
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Fig. 5. — The day —11 of the core-normal SN 2003du (solid line), from Stanishev et al. 
(2007), is compared with a synthetic spectrum (dashed line). 
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Fig. 6. — Three spectra of SN 1990N are compared: days — 14 and —8, from Leibundgut 
et al. (1991), and a previously unpublished spectrum obtained on day —12 by S. Benetti 
and M. Turatto. Vertical dashed lines refer to Si II A6355 blueshifted by 20,000 and 12,000 
km s" 1 and Ca II A8579 blueshifted by 25,000 km s -1 . The 6100 A absorption in the day -12 
spectrum appears to be a marginally resolved two-component (PV + HV) feature. 
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Fig. 7. — Spectra of one core-normal and four broad-line SN la of the one-week premax 
sample. Vertical dashed lines refer to Ca II A3945 and A8579 blueshifted by 25,000 km s~ x 
and Si II A6355 blueshifted by 20,000 km s -1 . 
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Fig. 8. — The day —7 spectrum of the extreme broad-line SN 1984A (solid line), from 
Wegner & McMahan (1987), is compared with a synthetic spectrum (dashed line). 
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Fig. 9. — The day —14 spectrum of the broad-line SN 2002bo (solid line), from Benetti et al. 
(2004), is compared with a synthetic spectrum (dashed line). 
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Fig. 10. — Spectra of one core-normal and three cools of the one-week premax sample. 
Vertical dashed lines refer to Ca II A3945 blueshifted by 20,000 km s" 1 and Si II A6355 
blueshifted by 15,000 km s _1 . 
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Fig. 11. — The day —5 spectrum of the cool SN 1999by (solid line), from Garnavich et al. 
(2004), is compared with a synthetic spectrum (dashed line). 
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Fig. 12. — Spectra of one core-normal and seven shallow-silicons of the one-week premax 
sample. Vertical dashed lines refer to Ca II A3945 and A8579 blueshifted by 25,000 km s~ x 
and Si II A6355 blueshifted by 15,000 km s -1 . 
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Fig. 13. — The day —9 spectrum of the shallow-silicon SN 2005cg (solid line), from Quimby 
et al. (2006), is compared with a synthetic spectrum (dashed line). 
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Fig. 14. — Spectra of one core-normal and three shallow-silicons of the early sample. Vertical 
dashed lines refer to Ca II A3945 and A8579 blueshifted by 25,000 km s' 1 and Si II A6355 
blueshifted by 20,000 km s _1 . 
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Fig. 15. — The day —12 spectrum of the shallow-silicon SN 1999aa (solid line), from Garavini 
et al. (2004), is compared with a synthetic spectrum (dashed line). 
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Table 1. The SN la Sample 



SN 


Epochs 


Galaxy 


References 




(days 









1984A BL 


7 




NGC 4419 


Wegner & McMahan (1987) 


1986G CL 


P. 




NGC 5128 


Cristiani et al. (1992) 


1989B CL 


7 




NGC 3627 


Wells et al. (1994) 


1990N CN 


— 8 - 


-14 


NGC 4639 


Leibundgut et al. (1991) 


1991T SS 


—7 - 


■1 1 

J. J. 


NGC 4527 


Phillips et al. (1992) 


1992A BL 


-6 




NGC 1380 


P. Challis, unpublished 


1994D CN 


-8, - 


-12 


NGC 4526 


Meikle et al. (1996) 


iyy toY bo 


-7 




T? Q P> K IP. C AH 


Li et al. (iyyy ) 


1998aq CN 


-8 




NGC 3982 


Branch et al. (2003) 


1998bu CN 


-6 




NGC 3368 


Hernandez et al. (2000) 


1999aa SS 


-7, - 


-12 


NGC 4469 


Garavini et al. (2004) 


1999ac SS 


-9, - 


■15 


NGC 2848 


Garavini et al. (2005) 


1999by CL 


-5 




NGC 2841 


Garnavich et al. (2004) 


1999ee SS 


-7 




IC 5179 


Hamuy et al. (2002) 


2001el CN 


-8 




NGC 1448 


Mattila et al. (2005) 


2002bo BL 


-6, - 


-14 


NGC 3190 


Benetti et al. (2004) 


2002er BL 


-7, - 


-11 


UGC 10743 


Kotak et al. (2006) 


2003cg CN 


-7 




NGC 3169 


Elias-Rosa et al. (2006) 


2003du CN 


-7, - 


-11 


NGC 1921 


Anupama et al. (2005), Stanishev et al. (2007) 


2005cg SS 


-9 




NGC 9290 


Quimby et al. (2006) 


2005hk SS 


-5 




UGC 272 


Chornock et al. (2006) 



Table 2. Fitting Parameters for Core-Normal SNe la of the One- Week Premax Sample 



rQ rnptpv 

-L CL1 Cl l 11C L CI 


SN 1990N 

Ul ' li/ Jul 1 ! 


SN 1994D 


SN 1998an 


SN lQQSbu 

Ul i U JU \ J LI 


SN 2001el 

Ul i ^UUILI 


SN 2003ce- 


SN 2003di 


v phot (km s^ 1 ) 


12,000 


13,000 


13,000 


11,000 


15,000 


13,000 


13,000 


v c (km s _1 ) 


1000 


1000 


1000 


1000 


2000 


1000 


1000 


r(C II) 




0.8 


0.4/14 










r(0 I) 




0.5/14 




0.6/12 




0.8 


0.3 


r(Mg II) 


1.2 


0.5 


0.4 


1.5 


0.5 


0.7 


0.5 


r(Si II) 


4 


5 


5 


8 


2 


8 


5 


r(HV Si II) 










0.5/20 


8 


5 


r(Si III) 


1.1 


0.6 


1.8 


1.8 


0.2 


1.7 


2.2 


r(S II) 


1.7 


1.5 


1.5 


2 


0.4 


3 


1.8 


r(Ca II) 


5(2) 


10 


4 


25 




25 


6(2) 


r(HV Ca II) 


1.8(4)/20 


2.5(3)/23 




1.6(2)/20 


20(8)/23[34] 


5/22 


5(3)/21 


r(HV Fe II) 


0.2(3)/17 


0.25(2)/19 


0.2(2)/19 


0.2(2)/17 


1/20 


0.4/19 


0.4/19 


r(Fe III) 


1.2 


0.5 


0.7 


1.5 




1.5 


0.8 



a For each ion the optical depth, r, is the optical depth at the photosphere or detachment 
velocity of the ion's reference line (ordinarily the ion's strongest line in the optical spec- 
trum). When a value of v e other than the default value in row 2 is used, the value (in units 
of 1000 km s _1 ) is given in parentheses. Minimum and maximum velocities (in units of 
1000 km s -1 ) are preceded by a forward slash. Maximum velocities are in square brackets. 



Table 3. Fitting Parameters for Three Core-Normal and Two Broad-line SNe la of the 

Early Sample 



L dl dlilt; tt:l 




qw i qqzLD 




QTVT 9nfl9hn 


QTVT 9nfl9pr 


Wp ho t (km s _1 ) 


14,000 


14,000 


14,000 


20,000 


15,000 


v e (km s -1 ) 


2000 


2000 


1000 


2000 


2000 


r(C II) 


0.5 


0.6 


0.6 






r(0 I) 




0.5 


0.3/15 


0.2(3)/21 




r(Mg II) 






0.6 


1 


1 


r(Si II) 


1.6 


8 


8 


30[25] 


6 [20] 


r(HV Si II) 


0.7(3)/20 




0.5(2)/20 






r(Si III) 


0.3 


0.7 


2.5 


0.6 


0.2 


r(S II) 


0.5 


0.6 


1.8 


1[22] 


1.5[16] 


r(Ca II) 




200 


12(2) 


20(10) 


10 


r(HV Ca II) 


8(6)/24[30] 


200(6)/23 


12(5)/22 




100(6)/23 


r(HV Fe II) 


0.3(3)/20 


1.2(3)/20 


0.3(3)/19 


l(3)/22 


0.7(3)/20 


r(Fe III) 


0.5 




0.8 







a For column descriptions see the note to Table 2. 



Table 4. Fitting Parameters for Broad-Line SNe la of the One-Week Premax Sample 



M T TY1 (^"1" PT 1 




SN 1 QQ9A 






w phot (km s _i ) 


16,000 


14,000 


14,000 


16,000 


v c (km s _1 ) 


2000 


1000 


2000 


1000 


r(0 I) 




1.5/16 


0.15(3)/15 


0.1 


r(Mg II) 


3 


1.2(2) 


0.8 


1.7(2) 


r(Si II) 


10(4)/[27] 


45 


7 


35 


r(Si III) 


0.6 


0.5 


0.5 


0.7 


r(S II) 


2.3 


2.5 


1.7 


3.5 


r(Ca II) 


20(6) 


50 


120 


150 


r(HV Ca II) 


20(6)/21 


2(20)/24[28] 


5(3)/21 


8(20)/24[29] 


r(HV Fe II) 


1.7(3)/24 


0.8/19 


1(3)/19 


1.5/20 


r(Fe III) 


1.4 


0.5 


1 


0.4 



For column descriptions see the note to Table 2. 



5. Fitting Parameters for Cool SNe la of the One- Week Premax Sample 





SN 1989B 


SN 1986G 


SN 1999by 


^phot l Km S ) 


1 o nnn 


1 o nnn 

1Z,UUU 


1 1 nnn 

11,UUU 


v c (km s 1 ) 


i nnn 

1UUU 


1 nnn 

1UUU 


1 nnn 

1UUU 




1 / ro 




A 


T-fMrr H 
T^lVlg ]_) 




9 ^ 
z.o 


l.o 


T{Mg LL) 


1 / ro 




o 
_ 




on 


An ( 9^ /h Ql 


fin /ri 4i 


r(Si III) 


0.7 






r(S II) 


2 


2 


1 


r(Ca II) 


10(2) 




40 


r(HV Ca II) 


l(2)/20 






r(Sc II) 






1.5 


r(Ti II) 






1.5 


r(HV Fe II) 


1/17 






r(Fe III) 


1.5 







a For column descriptions see the note to Table 2. 



Table 6. Fitting Parameters for Shallow-Silicon SNe la of the One-Week Premax Sample 







SN 1991T 


SN 1 997hr 


SN200 c ihk 


«N 1 QQQaa 


«N 1 QQQar 


SN 200 c irP' 


Vphot (km s x ) 


13,000 


10,000 


12,000 


5000 


11,000 


13,000 


15,000 


v c (km s _1 ) 


2000 


1000 


1000 


2000 


2000 


2000 


1000 


r(C III) 




0.3 




0.3 








r(0 I) 


0.1 


0.1 




0.08/6 


0.1/12 


0.1/14 


0.5/17 


r(Mg II) 


0.4(3) 


0.3 




0.35/6 


0.3(3) 


0.4(3)/14 


0.3(4) 


r(Si II) 


0.9 


0.2 


0.2 


0.2 


0.6 


1.8 


4 


r(HV Si II) 


0.4/19 












0.45(2)/21 


r(Si III) 


0.5 


1.3 


0.6 


0.3 


0.5 


0.9 


0.8 


rfS II) 


0.4 






0.25 


0.5 


0.4 


1 


r(Ca II) 


5 


1 


0.3 


0.8 


1.5(3) 


4(3) 


15 


r(HV Ca II) 


4(9)/20[28] 


0.6(3)/18 


0.2(2)/18 




1.3(9)/21[28] 


1.5(6)/19[23] 


5(20)/23[32] 


r(Ti II) 








0.04 








r(HV Fe II) 


0.5/22 








0.2/21 


0.3/19 




r(Fe III) 


0.6 


0.7(2) 


0.8(2) 


0.7 


0.6 


0.5 


0.8 


r(Co II) 








0.2 








r(Ni II) 








0.1 








r Ni III) 




0.1(2) 


0.15 


0.1 









a For column descriptions see the note to Table 2. 



Table 7. Fitting Parameters for Shallow-Silicon SNe la of the Early Sample 



± dl CllllC IjCI 


QM IQQQaa 


QN 1 QQQar 


kJl N ± U U ± ± 


v phot (km s -1 ) 


16,000 


17,000 


14,000 


v e (km s -1 ) 


1000 


1000 


1000 


r(C II) 




0.2/19 




r(C III) 


0.7 




0.5 


r(0 I) 


0.3/17 


1 




r(Si II) 


0.4/17 


5 


0.2 


r(Si III) 


0.6 


2.5 


2 


r(S II) 


0.5 


1.2 


0.5 


r(Ca II) 


1 


4(2) 


1 


r(HV Ca II) 


7/25 


1.5(5)/20 




r(HV Fe II) 


0.5/22 






r(Fe III) 


2 


1.5 


2 



a For column descriptions see the note to Table 2. 



